Close limit evolution of Kerr-Schild type initial data for binary black holes 
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We evolve the binary black hole initial data family proposed by Bishop et al. in the limit in 
which the black holes are close to each other. We present an exact solution of the linearized initial 
value problem based on their proposal and make use of a recently introduced generalized formalism 
for studying perturbations of Schwarzschild black holes in arbitrary coordinates to perform the 
, evolution. We clarify the meaning of the free parameters of the initial data family through the 

results for the radiated energy and waveforms from the black hole collision. 
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I. INTRODUCTION 



Collisions of binary black holes are expected to be one of the primary sources of gravitational radiation to be 
detected by interferometric gravitational wave detectors. Given the non-symmetric, time-dependent nature of the 
problem, the only realistic hope of modeling a collision is via numerical simulations. Unfortunately, given computer 
limitations, it is expected that in the near future most evolutions will have to start with the black holes quite close 
to each other. This brings to the forefront the problem of specifying initial data for the binary black hole collisions. 
Ideally, one would like to have initial data representing "astrophysically relevant" situations, that is, resembling the 
situation the two black holes would be in when they are in a realistic collision at the given separation. Unfortunately, 
_^ providing realistic data is tantamount to solving the evolution problem. Since this cannot be done, one is left with 
generating families of initial data based on mathematical or computational convenience. Initial attempts to provide 
initial data concentrated on solutions to the initial value problem that were conformally flat jjj. Conformally flat 
Q-( spatial metrics simplify considerably the constraint equations but suffer from drawbacks, most notably the inability 
to incorporate Kerr black holes, which are not known to have conformally flat sections (for a perturbative proof of 
non-existence, see |^]). Recently, attention has been drawn to the construction of initial data based on the Kerr-Schild 
form of the Schwarzschild and Kerr metrics. These constructions have several attractive properties: the slices are 
horizon-penetrating, which makes them suitable for the application of the "excision" technique for evolving black 
holes, and they can naturally incorporate boosted and spinning black holes ||. 

There have been two different proposals to use Kerr-Schild coordinates for binary black holes. In the proposal of 
Huq, Matzner and Shoemaker Q , two black holes individually in Kerr-Schild form were superposed. In the proposal 
of Bishop et al. Q the superposition was carried out in a way that the resulting superposed metric was in Kerr-Schild 
form. This latter proposal has the property that in the "close limit" in which the separation of the holes is small, the 
metric is given by a distorted Kerr-Schild black hole. 

In this paper we will consider the evolution of the Bishop et al. family of initial data in the close limit, by treating 
the space-time as a small perturbation of a non-rotating Kerr-Schild black hole. The evolution will be carried out using 
a recently introduced perturbative formalism that allows to evolve Schwarzschild black holes in arbitrary spherically 
symmetric coordinates . We will present an explicit solution to the initial value problem posed by Bishop et al. in 
the close limit and use it to compute the radiated energy in the black hole collision. This in turn will help use clarify 
the meaning of free parameters that appear in these families of data. 

II. KERR-SCHILD INITIAL DATA 

In their paper, Bishop et al. j|| assume that at the initial slice, the three metric and the extrinsic curvature are of 
Kerr-Schild (KS) type. The KS space-time metric is defined by 
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where fc M is null. The "background" metric i]^ is taken to be the Minkowski metric with coordinates (t,x) = (t,x l ) 
such that rjtt = — 1, f]u = and rjij — 5{j. The null vector fc M satisfies k t = —I and k l ki = 1, where k 1 = 5 13 kj. As an 
example, for the Schwarzschild geometry, 

M 

V = , hdx 1 = -dr. 

r 

For a KS metric, the three metric and extrinsic curvature with respect to a slice t — const, are 

fjij = 5. tj - 2Vkikj, (1) 
Ka = ~d t (Vkikj) + 2a [Vk s V s (Vkikj) - V (i (Vk 3) )] , (2) 

where a = (1 — 2V)~ 1 I 2 is the lapse and where V refers to the flat metric Sij . Bishop et aVs solution procedure 
consists in inserting (jl],||) into the constraint equations and to solve the resulting equations for V, V — dtV and 
ki = d t k il where fcj is assumed to be given. (k{ has only two independent components since k l ki = 1.) Below, we 
review the discussion of these equation where a particular ansatz is made for fc, representing two nearby non-rotating 
and non-spinning black holes. 



A. Two black hole data 



A single Schwarzschild black hole can be represented as 

h = ^ , |V0| 2 = ^ • v j( f>, 

with <f> = For two black holes, Bishop et al. make the ansatz 

M x M 2 
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where x_ x 2 denote the position of the black hole 1,2 which has mass Mi, 2 . If the black holes are located at = 
<Zi(0, 0, 1), with cii > > a,2, 4> may be expanded in a sum over multipoles: 

^ M ia [ + M 2 aj 

n r i = 2^ ^+1 p ^ cos ( 3 ) 

where Pg denote standard Legendre polynomials and where (r, 1?, tp) are polar coordinates for x. It is important to 
note that the expansion (Q) is only valid for r > max{ai, — a 2 }. 
Defining the separation parameter 

<2i — a 2 

e - 
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where M = Mi + M 2 is the total mass, and imposing the center of mass condition Midi + M 2 a 2 — 0, the close limit 
of (0) becomes 

M 9 MMiM 2 , , , 4x 

</»=—+ e 2 -^P 2 (costf) + 0(s 3 /r 4 ). 

As a result, to first order in e 2 , ki is given by 

k r = l, k A =e 2 ^^V A P 2 , (4) 
r 

where here and in the following, A = cp. The remaining amplitudes are expanded according to 

V = -— + e 2 v{r)P 2 , V = e 2 v(r)P 2 , k A = e 2 k{r)V A P 2 . (5) 
r 
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(In Bishop et ai.'s notation: v = vt and k — — fcy.) Introducing this into the constraint equations, and keeping only 
terms of the order e 2 , one obtains the equations 

3M / 2M\ • 3 6MM!M 2/ 
°=-- + — + — J fc --« ^(r-M), (6) 

, 4 6M 2 - 6MMiM 2/ 
=-v'--v + —k - (r - M), (7) 

2M- %MM X M 2 

= -MA;' + v + — k + . (8) 

Here, a prime denotes differentiation with respect to r. The system (0,|J) can be re-expressed as a single second order 
equation. Introducing the dimensionless quantities x = r/M and /i = M1M2/M , this equation reads 

= -« BX --« a , + - T t; + -5(3a: + 2). (9) 

Once we have solved this equation, the remaining amplitudes k and v are obtained from (|^) and (^), respectively. 



B. The solutions of equation 

A particular solution of ([)]) is given by 

2/j, ( 2 3 3 
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In order to find the solutions of the homogeneous equation, one performs the transformations x = 24/ z 2 , v(x) = z u(z), 
which yields the Bessel differential equation 

= z 2 u zz + zu z — (16 + z 2 )u. 

The solutions are a linear combination of the Bessel functions J&(iz) and Y&(iz). While J±(iz) behaves as z 4 for small 
\z\, Y±(iz) has the expansion Mj 

Qfi / 7 2 7 4 7 6 

Yjiz) = -—z- 4 — + - ^ + O(z 8 logz) 
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= -^-x 2 (l-- + \-\+0(x-Hogx)\ (10) 
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near z = 0. Thus, the general solution to (^|) is v(x) = fiv(x), with 



= -|(i_| + J + |. + ^r 4(l y24^)) + g j 4(jv ^). (ii) 

While u is regular at x = 00 for any values of the constants C\, C2, equation (Q) shows that in order for to be 
regular at x = 00, it is necessary that u decays at least as fast as x~ 2 . Comparing the expansion ( jUj ) with ( O ) 
one sees that by choosing Ci = 67r, one can get rid of all terms which decay slower than x~ 3 . By looking at gauge- 
invariant expressions, we will show later that this choice is indeed necessary in order to get an asymptotic flat solution. 
Therefore, C\ is fixed by physical means. The role of the constant C2 is discussed below. 



C. The Zerilli amplitudes 

In Ref. H , we have recently derived a gauge-invariant generalization of the Zerilli equation which allows to study 
perturbations on a Schwarzschild background written in any spherically symmetric coordinates. In Appendix A of 
this paper we have written the perturbed metric in terms of the generalized Zerilli function ip, that one obtains using 
the formalism of |^| for the case at hand: I = 2, even parity perturbations of a KS background. This metric is a 
solution of Einstein's vacuum equations provided the Zerilli function xjj satisfies 
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where ip = "0( r > T : = V-^) an( ^ now ^ = "0' = ^x^- 

In order to evolve the KS initial data, we have to relate the amplitudes v, v and k to the scalars ip and ip (introduced 
in ||) which satisfy the Zerilli equation (|l^). Using the expansions (Q,^) in the expressions |l]J^), it is straightforward 
to calculate ip and i/> using the formulae given in §. The result is 



ip = M/i 
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where we have also used (Q) in order to eliminate k. Since ip is gauge-invariant, it is clear that the free constants C\ 
and Ci appearing in v(x) cannot represent a gauge freedom. In order to have an asymptotic flat solution, v(x) must 
vanish at infinity. As discussed in the previous subsection, this fixes the value of the constant C\. In this case, ip 
and ip fall off like x~ 2 at infinity. The constant C2 is still free and will determine different sets of initial data as one 
chooses its value. The radiation content, as we will see, depends on C2. The constant is therefore clearly associated 
with the "spurious radiation" that the initial data contains with respect to "astrophysically relevant" initial data. 
One could probably determine this content by evolving the initial data set backwards in time. This calculation would 
be possible (at least for a limited amount of time) within the confines of the close approximation if the black holes 
are initially very close. One could therefore follow the space-time backwards for a short time and see if incoming 
radiation is present at a finite distance of the holes. We have not performed such a study, but it is feasible (we thank 
Jeff Winicour for bringing this to our attention) . 

It should be noticed that the initial data for the Zerilli function diverges in the limit x — > for all values of C2, so 
one cannot single out a preferred value of this constant by demanding the initial data to be finite in this limit (even 
though, as already mentioned, the multipole expansion is, in any case, valid only for r > max{ai, — &2}). 



D. The linearized apparent horizon equation 



Bishop et al. have argued that the position of the apparent horizon is related to the constant C2. Here wc perform 
a linearized analysis of the position of the horizon, to clarify the meaning of their finding. 

Given initial data Qij, Kij on a space- like slice S, the location of an apparent horizon (AH) can be determined by 
the equation 

V, t s l - K ij s i s j + K = 0, (14) 

where s l is the unit outward normal to the AH. If the AH is given by f(x) — for some function / on E, we have 
Si = AVi/, where A" 2 = V l /V;/. For KS initial data of the form ([l]J^), it was assumed in |5] that the AH coincides 
with a surface S which is orthogonal to ki . In this case, = —ki/ot, and one can check 
compatible with the result in ||, i.e. S is an AH if V — —1/2 on S. 
For spherically symmetric initial data, 



c that equation (|i4|) is 



gijdx 1 dx° = ^{rydr 1 + r^dQ, 



Kijdx l dx j = p(r)^(r) 2 dr 2 + q(r)r 2 dVt 2 . 



we must have s r = 7, sa = 0, and the AH equation (|14| ) yields q = —l/{r^f). It is not difficult to show (either by 
using the KS form of the Schwarzschild metric or more generally by integrating the constraint equations) that this is 
equivalent to r — 2M, where M is the ADM mass. We now want to linearize the AH equation around a spherically 
symmetric background and to find the deviation of the AH from r = 2M . In the linear regime, we expect that the 
location of the AH can be described by the image of the circle |x| = 2M under a map of the form 

x 1— > x — e 2 D(x) = . 

r 

The deviation function D(x) is related to the function f(x) as follows: 
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= / (x- e 2 D{x) = ) = + e 2 [sf(x) - 9 r / (0) fe)D(x) 



where f^°'(x) is a function describing the AH to zeroth order (for example f^°'(x) — r — 2M). Using also A = 7(9 r /)~ 1 , 
we obtain D(x) = XSf/j. 

On the other hand, for linear perturbations around a spherically symmetric background, the normalization of s l 



yields Ss r = 7/1/2 while equation (14) gives 



r 2 



= V A (jS SA -q A ) -rh+—k' + 7 r 2 V k , (15) 



where h, k, q A and Vk are defined by 

h = ^ 2 5g rr , qs = 5g rB , k = g AB Sg AB , V k = S(g AB K AB ). 
In terms of the function /, s A — 5(Xd A f) = -fd A D, and the linearized AH equation finally becomes 

r 2 

7 2 AL> = V A q A + rh- —k' - jr 2 V k . (16) 

Performing a multipolar decomposition of D(x), this equation becomes a set of algebraic equations for D. Evaluating 
for the KS data proposed in Section II A, we find that 

D(x) = a 2 I — (3r + 2M)t-(r) J P 2 (cos^). (17) 

So we see that the position of the apparent horizon is given by the image of the circle x — 2 under the map 

x i-> x — s 2 D(x) — 

X 

where the deviation function -D(x) can be expressed algebraically in terms of the perturbed three metric and extrinsic 
curvature. For a KS metric and x = 2, 

B(£) = |Q«(2)-l)ft(cos1?). 

It is now clear that the deviation function depends on the value of the constant C2 . In particular, we can choose C2 
such that D{x = 2) vanishes. 

A way to see that the meaning of the constant C2 is not just a choice in the position of the apparent horizon is to 
notice that once one has fixed the values of M% and M2 , the KS form of the metric completely fixes the coordinates 
on the KS slice (at least to linear order). Indeed, a gauge mode must satisfy the constraint equations (|^,|^,||). On the 
other hand, the general solution to these equations is completely determined by v(x), which cannot contain gauge 
modes since it is related to the gauge- invariant amplitude tp according to (p"3|). Therefore, we have two possibilities to 
fix the coordinate location of the apparent horizon (at x = 2, say): The first possibility is to perform an infinitesimal 
coordinate transformation such that the apparent horizon appears unperturbed relative to the Schwarzschild horizon. 
In this case, the initial data will not have KS form anymore. The second possibility is to adjust the constant C2 
such that the apparent horizon is at the location we desire. Clearly, these two methods are different since the former 
corresponds to a gauge transformation, while the latter corresponds to a true physical change in the initial data, as 
we anticipated before. So indeed the constant C2 is related to the position of the apparent horizon as was noticed by 
Bishop et at, but through a genuine change (not just a gauge change) in the initial data. 

Figures 1 and 2 show the initial data for the Zerilli function and its time derivative for different values of C2. In 
the next section we shall analyze the dependence of the total radiated energy and waveforms on C 2 . 



III. EVOLUTION 

An expression for the radiated energy in terms of gauge- invariant quantities is given in H. Since here we have 
expanded all perturbations with respect to the Legendre polynomial P2(cos$) = Y 20 ({)), this energy expression 

becomes 
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FIG. 1. The initial value of the Zerilli function and its time derivative, 
choice of scale of the figures. 



C2 = is nonvanishing, as may appear due to the 



dE 6 



du 



with tp being evaluated in the radiative zone. 

We have written a code that solves our generalized Regge- Wheeler and Zerilli equations. As a consistency check, we 
have evolved the close limit of some maximally sliced initial data (which can be seen as perturbations of Schwarzschild 
in usual coordinates), being able to reproduce previous values for the total radiated energy (e.g. Misner's initial data 
||, or boosted black holes 0). 

The code is a standard second order dissipative, finite differencing, one. In the case of a KS background we perform 
excision, i.e. we place the inner boundary inside the black hole, and in that way avoid giving boundary conditions 
there. In figure 3 we show the Zerilli function, scaled by /i (i.e. ip/fi) versus time, extracted at r = 100AT for two 
different values of C2. From that plot one can notice, for example, the typical ringing frequency for Schwarzschild 
black holes. 




C2=0 

C2 = 100 



50 100 150 200 250 300 

t/M 

FIG. 2. The radiated waveforms at r = 100M for two different values of C2. As usual in close limit collisions, the waveform 
is dominated by the fundamental quasi-normal mode. 

Given the linearity of Zcrilli's equation, and the form of the KS initial data, the dependence of the Zerilli function 
on the parameters of the problem is 



^(i, r) = e 2 /x (j} a {t, r) + C 2 ^,(t, r)) 



where the functions ip a and tpb are dimensionless. Accordingly, the total radiated energy is 

A, ,2 / /-oo 



E = 



6eV 



i>a dt + Ci / tp b dt + 2C 2 \ i> a ipbdt 
Jo Jo 



Therefore one needs to perform only three runs to obtain the complete dependence of the radiated energy on the free 



G 



parameters. The result is 

E = e i Mn 2 (3.6 x 1(T 4 + 5.2 x lCr 7 C 2 2 - 2.2 x KT 5 C 2 ) 

The quantity -E/(e 4 M/i 2 ), as a (quadratic) function of C2, has a local minimum at C2 ~ 21, where it takes the value 
E/(e 4 Mfi 2 ) w 1.3 x 1CT 4 . This function is plotted in figure 3. 

Total radiated energy 
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FIG. 3. The radiated energy (in units of e 4 Mfj, 2 ) as a function of G%. As a rough comparison, if we consider equal mass black 
holes (fj, = 1), and take the "separation" e = 1 and identify it with the "separation in the conformal background geometry" 
for two Brill-Lindquist black holes (for a more physical picture, a conformal separation of less than 0.8M corresponds to a 
common apparent horizon for the Brill-Lindquist family), the latter would radiate around 10 -5 M, which is rougly similar to 
the radiation we get for the minimum value of Ci . 



IV. DISCUSSION 

We have evolved the initial value family of Bishop et al. in the limit in which the black holes are close to each other 
by treating the spacetime as a single distorted Kerr-Schild black hole and solving the linearized Einstein equations 
for the distortion. The evolution sheds further light on the role of the integration constants present in the family. 

An obvious question to ask would be "does this family contain more/less radiation than other families" (for instance 
the Misner data). Unfortunately, the family has explicit free parameters and therefore the comparison is highly 
dependent on the arbitrary values of these parameters. This should not be misinterpreted as a problem: it just 
highlights that the initial value problem for binary black holes inevitably contains ambiguities. Some proposals may 
resolve the ambiguities based on aesthetic criteria, but from a physical point of view that is not more satisfactory 
than simply picking values for the constants involved. These issues could be better understood if one evolved the 
systems backwards in time and tried to establish the amount of incoming radiation. The present results should be 
of interest in the calibration of numerical codes based on the Kerr-Schild coordinate system. Experiments with the 
Maya binary black hole code ]Tc| to compare results are currently under way. 
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APPENDIX A: PERTURBED METRIC FOR A KS BACKGROUND 



For the particular case in which the background metric is KS, the perturbed I — 2 even metric in the Regge- Wheeler 
gauge, in terms of the Zerilli function, is 



1 2 6 
g rr = H h — 

x M 



6(x + 2)(3 + 6x + 4x 2 + 4x 3 ) , 4(4x 2 + 9x + 3) • 



(3 + 2x) 2 x 3 



x 2 (3 + 2x) 



2(2x 3 + 2x 2 + 15a; + 6) , 



x 2 (3 + 2x) 



if)' - 2x0 



an 



12(3 + 6x + 4x 2 + Ax 3 ) , 4(2x 2 - 6x - 3) , , 



2 

x + M 
2(x + 2)(2x 2 - 3 -6a;) 



x 3 (3 + 2x) 2 



x 2 (3 + 2x) 



x 2 (3 + 2x) 



— 2x0' 



9eo = x +- 



12(x + 1 + x 2 ) 
3 + 2x 



■0 - 4x0) - 2x(x - 2)0' 



K 20 (tf) 



gee sin 2 



fftt = -1 



x M 



6(4 + x 2 )(3 + 6x + 4x 2 + 4x 3 ) 4(5x 2 + 15x + 6) • 

~^ + x 2 (3 + 2x)(x + 2y 



(3 + 2x) 2 x 3 (x + 2) 



2(2x 4 - 2x 3 - 5x 2 + 24x + 12) x - 2 8x • 

2/r, , 7TJTFZ i r>\ Y ~~ %X - ; o — o 



x 2 (3 + 2x)(x + 2) 



x + 2 



x + 2 



where 5 is a perturbative parameter. 
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